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Span load di stri out ions of sv/ept-'bacic, v/in^s have "been 
calculated. The neth.od used was to replace the wing v/ith 
a* "bound vortex at the quarter-chord lino and to calculate 
the downwash due to the s^^ston of ootmd and trailing" vor- 
tices to conforn a.t the three-quart cr~chord line to the 
slope *of the flat-plate wing surface. Hcstilts are ji^iven 
for constant-chord and 5:1 ta"OGrod iDlan forns, for sv;eop- 
hack angles of 0", 50 , ■ and 45 , a^nd for aspect ratios of 
3, 6, and 9. Sonc connonts on the stalling of swopt-*back 
wings are included, 

IlTTEODUCTIOil ■ 



Airplane desi^^ns are often proposed in v/hich stability 
and control are provided hy a v'ing with a largo angle of 
sv/ecpoack instead of "by the conventional fuselage and tail 
arrangenont , The calculation of the aorodynanic forces is 
noro difficiilt on such a wing than on the corresponding 
unsv/ept-hack wing "because of the nore conplicated nature of 
the vortex sjrsten, the flow of which is to ho calculated. 
The prohlen is essentially "that of tvio yawed airfoils (the 
two halves ox the swept-hack v/ing) subject not only to the 
effects of finite aspect ratio "but also to their nutual 
interference , 

Ivieghardt has shown in his work v/ith rectangular flat 
plates (reference 1) that good accuracy in the total lift 
force can he attained hy replo.cing the plate v/ith a "bound 
vortex at the quarter-chord line and naking the dovmvfash 
due to the "bound, vortex and its associated trailing vortex 
systei.i conforn, at the three-quarter-chord line, to the 
slope of the flat plate ♦ 



Such a procedure is rigorously correct for an unyawed, 
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rectangular flat plate of infinite span. The nethod also 
holds for a yawed, infinite-span fla.t plate "because, in 
this ca.se, the spanv;ise conponent of the free— strea.n ve- 
locity has* no effect ( if. potent ial flov/ is assuned) on the 
lift v/hereas the conponent nonial to the span has the sane 
effect cas on a.n unyav;ed airfoil. More precisely, let the 
airfoil 'be yav/ed through an angle A at zero angle of at- 
ta.ck and then rotated through an .angle a o-hout the line 
with respect to v/hich A was neasured. If the free— strean 
velocity V is resolved along the throe nutua,lly perpen- 
dicular directions - nar^iolj^, spanwis^e, chordwiso perpendic- 
ular to the span, n.nd nornal to the plane of the flat 
plate - the velocity conpononts arc, respectively, 
y sin. A cos a, V cos A cos a, and Y sin'a. Inasimch as 
the chord neasured norna.l to the span is c^ cos A| then 
fron t wo-dinensi onal thih-airfoil theory-, the circulation 
r is given "by 

r c rr ViCq cos A^ sin a (l) 

Also, fron t wo-dinensional wing theory, the center of 
pressure of the lift forces is at the o^uart er-chord line. 
If a hound vortex of circulation given' "by equation (l) is 
placed at the quarter-chord line, the nornal induced veloc- 
ity w at the thr ee-qua,rter--chord line is 

p TT V Co COS A sin a 

• w = — — = 

2Trr . « Co 

2tt — cos A 

or 

w = 7 sin a (2) 

which is exactly the anount . ne eded to."produce a resultant 
flov/ along the flat-plate equivalent of the. thin airfoil 
at the three-quarter-chord point . 

The extension of V^ieghardt * s nethod to flat plates of 
swept-hack plan forn has heen nade in the. present paper. 
A swept-hack hound vortex was placed at the quarter-chord, 
line (fig* l) and the dov;nwash due to the hound vortex and 
its trailing-vo:tvfe-e-'X sy sten \iras calculated at points on the 
three-quarter-chord liaie. .The "bound vortex strength was 
deternined Dy the condition that the resultant flow at the 
three-quarter-chord line he tangent to the flat plate. 
The span load distrihution thus * det ernined should , he ex- 
pected to he less valid near the. tips and the center of the 
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v/ing' thsji in between. In the cent or , .-the . effect • 6f one v/ing 
upon the other is strong and the hound vortices are curved 
lines rather than hrohon lines- v/horeas, at the tips, the 
presence of strong trailing vortices invalidates ■ thc^ argu- 
nent of the previous para.graphs, which was based on two- 
dinensional wing theory* : 

The results obtained -for the flat-plate wing surfaces 
nay be applied to wings having the usual type of airfoil 
section on the assunption that the finite thickness and 
canber of the wing do not appreciably displace the vortex 
systen fron the horizontal plane. lurtheruore, inasnuch 
as the section lift-curve slope aQ generally differs fron 

2-17, the airfoil chord to be used with the results of this 
paper is 

^0 • ' . 
X c *" — — c " 

2tt c.ct\ial 

This correction is based on the -t wo-dinensi onal- airfoil re- 
lation 

V ' 1 

2 ■ ^ 

in which ca^ occurs as a'product* The value of a^ v/ill 
depend on the angle of swee-pba.ck A a,s well as on the air- 
foil section. If there v;oro no vi scosity of f ect s nor 
trailing-vortcx systen, a^ would be deternined by the. 
section perpendicular to the span because the conponent of 
the air flow perpendicular to tht) span produces the lift. 
The chief cause of the reduction' in a^ fron its theoret- 
ical value as given by potential theory is, however, the 
viscosity phenonena (boundary layer, separation, etc.) that 
depend on the. actual path of the air particles over the 
airfoil. Because, for largo aspect ratio at least', the 
streanlines over the airfoil do not greatly depart fron 
the frce-strcan direction, ' the slope of the lift curve of 
the airfoil section in the free-strean direction should be 
used, • • • 

SYIOOLS' 



The * following synijols . are used: 
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"b vfi-ng span, noasurcd along q^uart or-chord lino fron 
tip t'o tip 

8 v;ing sonispan ("b/s) 

A sweepback angle, neasurod fron quart cr.-ciaord lino 

a angl'o of attack, noasured in free stroan direction. 

S. area of ' plan" f orn of v/lng 

c v/lng. chord neasurod in- f roe-strean* dire ct ion - 

Cq v/ing root chord noasured in froo-stroan direction 

-^A aspect ratio g -J 

^-L total lift 

a^ slope of section lift curve 

lift coefficient ( ^ 

p ^' o 

I lift per unit length along span 
r circulation 
Dj_ induced drag 



C-py. - inducod-^drag coefficient 

P 



'1 



p density of._air 

V frce-strcan air velocity . 

V/ induc-ed dov/nv/ash. velocity . 

Coordinates: 

T| coordinate of downv;ash point in direction of quarter 
chord line (positive to right) 

a coordinate of TDound vort ox . e-lonent in direction of 
quarter-chord line (positive to right) 

or center of pressure of span load distri"bution 
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J8 • angle corresponding, to.. rO". (cos 0 = 

\|/ angle corro spending to T) (cos \|/ = r|/.s) 

y perpendicular to free-stream direction in plane of 
v/ing (positive to riglit). 

Subscripts: >^ 

+ refers to right side, .of v/i-ng 

- refers to left side ox wing 

1 refers to hound vortex 

a refers to trailing vortex 

laETHOD 



The flat-plate -v/ing surface of swept-hack plan form 
is replaced v/ith a oound vortex at • the quarter-chord line 
(fig» 1). The hound vortex and an element of the trailing- 
vortex- sheet as-- ;V/ell- as the svmhols us.edj In- .th:e- following 
expre.ssi-on:s,. ar:e' shov/n .in figure 1.* By- the .Biot^-Savp^rt' 
law the • induced/ velocity w (positive downward) . cat. the 
point P- on the three-quarter-chord lino, due to the 
bound vortex r on the quart§.r-chord line, is .. ^ 




^^^1 = / " ~ + — / cia (3) 



which, -upon substitution of the'^ geom-e'trical relations 
cos e+ = - - ^ sin e+ = ^ 



- Tl + B tan A 3 
cos e_ = = = sin 9_ = — 



(4) 



become s 
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Amij^ = - ^ J* r (e + ) sin 6+. d6+- 



'+0 



- ^ / r(e ) sin e_ de_ 



(5) 



The corresponding expression for the trailing vortex sheet 
is 



s 
n 



3 



1 COS 



(e, - A . I)] dr 



0"^__cos Aj- y 



0 r 



1 + cos (e_ + A - 



dr 



a cos A - y 



(6) 



-s 



The total induced velocity is the sub of equations (s) and 
(6); this sum "becomes, after equation (5) is integrated "by 
parts and it is noted that F is zero at the" wing tips, 



4„,, = + • J- ^ - 



cos 



e+ dP 



1 



cos e_ dP - 



[1 + sin (e+ ^ A)3 
cos^ - y 



dT 



0 



r 

-J 



-s 



1 + sin (6^ + A) 
a cos A - y 



dr 



(7) 



The variables of integration a, 6, and T may lie ex- 
pressed in terms of a single independent varia*ble 0,. 
varying continuously from 0 to rr across the span from 
right to left. Thus hy the definition 
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^•^ = s cos 0 for -s <.a^ < 0 

■■■ .: :- . (8) 

= s COS 0 for 0 < (J^ < s 
and the geonietric relation. ' , 

tan e. = — tan 6' = = — (9) 

I 

together v/ith the usual Fourier series expression for T 

00 

r = 4itVCo sin a Z a^,^.^ sin (2k + l) 0 (lO) 

k=o 

the follov;ing expression may he derived hy substitution 
into equa-tion (7): . . 



^ . . = E (2k+l) apv+i 
Vcq sm a ■~-"^-** 



2k+l 



TT/2 ' TT/2 

+ /' COS (2k4-l) '0 d0 _(A .^Q^ 0 - A4>) cos (2k4>l)p 

s COS A COS 0 y ' ^+ J y + (s cos 0 - JLj.)^ 



+ B4. cos A 



TT/2 . ' 

■ . cos (2ic+l) 0-d0 



(s cos A cos 0 - j) J "^ty + (s COS 0 - kj^y 



Tr/2 

(s COS A - Af) COS (2k+l) 0 d0 
- sin h J : ^ — 

0 



(s cos A cos 0^- y) J '£>j^ + (s cos - A4.)' 



TT ' . . _ • . _ Tt 



+ / cos (2k+l) 0 ^0 ^-J^ P {b cSs 0 A^) cos (2k+l) 0 

J s cos A co^ 0 ^ y 3^ J / 4- (s cos 0 - A )^ 

TT/a «- TT/2 V - r ^/ 



TT 



+ B COS A 
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cos (2k+l) 0 d0 



n/2 



3„ + (s cos 0 - A^)" 



TT 

+ sin A ^ (s cos 0 - A,) cos (2k+l) 0 d0 



TT/2 



(s cos A cos 0 - y) y +' (s cos. 0 - A_)^^ 



(11) 



The first and the fifth integrals .are integrable "by means 
of 



^ CQS 31 0 TT 



sin n \|/ 



cos 0 -^c•os \|/' sin \|/ 



(12) 



The singularities in the integrands of the third, the 
fourth, the seventh, and the eighth integrals can T^e removed 
"by the identity ' 



cos n 0 



(.os0-5)/l.[i-(cos0-^)" 



V (A^ -Ti)^ 



COS n /) 



cos 0 - -Q. 

^ s 



A 

COS n0 (^cos 0 + ^-2 



yi.r|-(c6s0-.^) 



Ti 

.1^ 



/3/*(A+-'n) +/i + 



1 



'7J 



together v;ith the elomchta,ry integration 



(13) 



TT/8 

rr ^ r : cos n0 d0 

J cos 0 - 'cos \l/ 



(14) 
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The limits 



r 

s / can 
TT/a 



"be* changed to /' hy the substitution 

^0 



0 = TT - 0 . VJith these 'cha^nges and the expression of all 
distances in terias of the root chord c^., equation (ll) 
for the dov/nwash at P "becomes 



CP 

= E (2k+l) a, 

Y sin a k=o 



■2k+l< 



(-1) 



k / 



2k + 1 



TT/a 



^ TT sin (2k+l)\i/ ^ s 
s cos A sin \i/ 



cos (2k+l) 0 cos 0 



I- 



d0 



TT/2 



1 /Af 



r cos (2k.l) 0 ^ , ^^ -^ 



TT/3 



K. 



2k+l 



-(i)7 



cos 



(2k+l) 0 (cos 0 - B > 2 I) 
P (C+J) 



TT/2 



s / COS 0 cos (2k+i) 0 



d.0 



TT/a 



TT/a 

T)) tanA\ r cos (2k+l) 0 

(cos0^-n> 



where 



(15) 
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- (cos 0 - -| 



1 + 



c = 



3. 



8 = A . - Tl. 



(16) 



COS \|; 

and where 



n 

s 







K3 = 


4I3 - .3Ii 


Kg = 


I6I5 ^ 2OI3 + 5Ii 


= 


„64Iv - II2I5 + 56I3 
etc , 


lo = 


• 1 , loge ^ ^ 
sin \i/ cos \1/ 


Ii = 


— + cos \^ Iq 


^& = 


1 + cos \J/ I J 


I3 = 


J + cos \if I2 


I4 = 


— + cos \^ I3 

3 . 


I5 = 


3lt + cos \j/ I4 
16 ■ 


Is - = 


. + cos Aj; Is ■ 


I7 .= 


|5 + cos \|/ le 
32 



(17) 



/ (18) 



etc . 
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The integrals in eq.uation (15) can te expressed in 
terms of' elliptic functions tut the' formulas are too long 
for practical evaluation. The integrals were therefore 
evaluated numerically hy Simpson's rule and the downwash 
was thus computed as a linear function of a^, a3, as, 

and a^ at four points on the three-quarter-chord line. 
These points were located at \j/ = cos"^ ^ = 20*^, 40°, 60°, 

and 80° corresponding to — = 0,940, 0,766, 0.500, and 

s 

0.174» The condition that the resultant flow "be along the 
flat plate at these locations (and the symmetrically .sit- 
uated spanv;ise positions) was fulfilled hy setting 
V/ = V sin a in equation (15). F.our simultaneous equations 
were thus ohtained and solved for a^ , a3., a^, and a^, 
v/hich hy equation (lO) gave the spanwise rdi st r ihut i ons of 
circulation. The convergence, of .the . Four ier series- thus 
ohtained for the circulation r was rapid. (See tahle I.) 



RESULTS 



Span load distrihutions Were calculated, for constant- 
chord and 5:1 tapered plan forms with ''3'0° angle of sweep- 
back and for constant chord plan form v/ith 45° angle of 
sv/eephack. The aspect ratios calculated v;ere 3, 6, and 9, 

The Fourier series coefficients for the circulation, 
the spa;nwise locations of the center of pressure, and. the 
induced-drag coefficients are given in tahle I» The 'cen- 
ter of pressure of the span load distribution is given by 

s tt/s 

^ CO 

apvr cos A da / cos 0 sin 0 S ^k+1 (2k+l) 0 d0 




s tt/s 
s / pVPcosA da . / sin 0 L ^2k+l (2k+l) 0 'd0 



o 0 



4 / 3 a, 1 a^ 1 \ 

- — (l+--^-- — + T^— • .) (19) 

Stt \ 5 a^ 7 a^ 15 a^ / 

The induced drag was obtained by projecting the lift- 
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ing. ele-nents -into .one line perpendicular-' t.p trhe f r ee-st r ean 
direction (Hunk^s stagger theoren, reference -2) and olstadn- 
ing the induce.d- drag of the .re sul ting - llf t ing-line distr.i- 
"butiqn.. Thus, for 'the distribution . ; - 

.. r = 417. V:.Cq .sin ..a 2 ^2k+l ^^^^ (2k+l.) 0 

• ' ' 0 . [ . • . , 

regarded as a lifting-line, distribution, the induced drag is 

" s c 0 s A ' ' ' 

1 / . V cos A 

0 

Tr/2 

Z (2k+l) a2k+i sin(2k+l)0 £ agn+l sin(2n+l)0 d0 



0 



o 



With the relation "betv/een Cj^ and a^ derived in .equation 
(21) • ^ 

^ = Z ..(2k + l)/^2^) • '(20) 

ttA 0 V a^^ / 

which holds for both the constant-chord and the tapeired 
plan forns, 

The. total lift forces, are obtained fron a£ as fpl- 

lov/s: 

Constant chord: * • ^ 

T ■ P V cos A / r ds 
Cy = L_ = 1 :z| : — = 2tt^ a ai 

t P S p S 

5:ltaper;, 

Cl = ; — = -3- tt2 a. ai 



cos a 



(21) 



The variation of the total lift forces with aspect ratio 
• is given in figure 2. - The curve..- for zero- angle- of sweepliack 
is practically identical with that given Tsy V^ieghardt (ref- 
erence 1, fig. 4) "but is extended to as aspect ratio of 9, 
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Included for comparison -are the curves for the correspond- 
ing unsv/ept-ha'ck plan forms olDtained from lifting-line' 
theory, a curve due to Krienes (reference 3) for an ellip- 
tical flat plate calculated on the "basis of the accelera- 
tion potential, and a result from experimental v/ork (ref- 
erence 4) for a 2:1 tapered v/ing of 30^ sv/eephack. 
There is a surprisingly large difference het'w'een the total 
lift" forces predicted "by the ordinary, lifting-line theory 
and those predicted "by the modified theory used in this 
report. The experiments of V/inter (reference 5, fig. 36) 
on rectangular flat plates favor the modified theory. 

For a low aspect ratio and an ar'bitrary plan form, 
equation (24)' for the limiting span load distribution 
(presented later in this paper) yields 

. = A (22a) 

A— ^-0 27701 4 

while for large aspect iratiq and ar'bitrary taper, equation 
(l) ■ gives 

A^^"L ^ = A- . (22h) 

The reduction of lift "by cos A, indicated "by the preced- 
ing expression, may "be regarded as a first approximation 
to the lower lift force of a swept-back v/ing» Figure 2 in- 
dicates that the approximation "becomes loss valid for lower 
aspect ratios, higher angles of sweepback, and higher 
amounts of taper • v. * 

A comparison of.-the curves for the constant-chord, 
30.° and 45 sweepback plan forms,' and the .r.o-ct angular plan 
form (by Wieghardt * s method) shovfs that the total lift- . 
curve slopes are reduced loj about. 7 -and 19 percent for. the 
30° and the 45° cases, respectively, a,t an aspect ratio of 
6. If the total lift of the rectangular plan form ob- 
tained by lifting-line "theory is., used as. a. basis of com- 
parison, the corresponding ..tot a.1 lift reductions .are 14 
and 25-,percent, 

Figure 2 also shows that the ..rcduct i on in lift- of the 
swept-back taporod. plan form^ is much less than that^ of the 
swopt-back constant-chord,.. plan .form. 

The span load di stributions. arc s'hov/n in figure 3. 
The ordinates are given by 
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(23) 



It - can be shown t;hat , as the aspect- .ratio- approachos..'ze;ro, 
equation (ll) -yields-,- for -an ■ arh'i trary plan. form, 

^ •■ ' • ■ * ' .. 

lini . • ^ L^^ _ = '± ' stn. 0' • ■ (-24) 

* A— I p T s Cl ^ ' ■ ■ . 

as a limiting form. 

Figure's sho.w.s that the wings with the swopt-hack 
constant-.chord' plan forms have higher tip loadings than 
•the corresponding unsv;ept-back wings." This result implies 
greater tip-stalling" tendencies for sv/cpt-hack v/ings. 
Furthermore, the root sections tend, with higher aspect 
ratios, to "bo less loaded than the sections halfway botweeii 
root and tip. That this result must- be valid is evident 
from the fact that a yawed, infinite-aspect-ratio, constant^ 
chord airfoil has a uniform span lovad distrilDution. If one* 
half of the airfoil is hent back to form a sv/ept-back wing 
and if the hound vortex strength remains the same, points 
near the root section, which originally had the proper 
downwash to conform to the slope of. the surface, would now 
have too much down\^rash; hence the bound vortici ty must be 
reduced- at and near the root sect ion • 

Figure 3(a) indicates that lifting-line theory pre- 
dicts higher tip lo'ading than is given hy the modified 
theory used herein. llov/, the trailing vortex system in- 
duces ..less dov/nwash at the quarter-chord line than at the 
three-.q.uartor-chord line^ Since lifting-line theory de- 
ducts the dov/nwash angle at the quarter-chord line from 
the geometric angle of attack," the lift near the tips 
therefore appears higher than that given by the modified 
thoo ry » 

The spanv/ise location of the beginning of the stall 
on a sv/cpt-back wing and the corresponding angle of attack 
are, hov/evcr, primarily dependent on the spanwise flow of • 
the boundarjr layer on the 'suction surface. This spanwise 
flov/ is due to the surface pressure .distribution*. On a 
swept-back wing' th"© surface pressure gradients si-zeep the 
slov/er moving air of' the "boundary I?:yer "toward the' tip. 
Some tuft observations on con.st ant-chord swept-back and 
swept-f orward wings indicated 'that ,, in addition,.^ the tip 
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vortex, which was coining off the top of the wing, provided 
an added local inducement for the boundary layer to sweep 
toward the tip. Near the tip, the thicker boundary layer 
(before it has been excessively influenced or removed by 
the favorable pressure gradient of the' tip vortex) will 
therefore stall the wing first in that region. On a swept- 
forward wing the surface pressure gradients sweep the 
boundary layer toward the center but near the tip vortex 
the flow is still outward because the tip v.orte'x comes off 
the top. of the wing. 

The span loadings of the tapered plan forms do not 
differ greatly either in magnitude or distribution from 
those that lifting-line theory predicts for the correspond- 
ing unswept-back wings* The total lift forces on the 
tapered plan forms^being higher than on the c ons tant— chord 
plan forms, stalling should occur at a lower angle of at- 
tack for a tapered wing than for the corresponding constant- 
chord wing. If two— dimens ional chordwise load, distribution 
on each spanwise location of both the tapered and the 
constant— chord plan forms is assumed, geometry indicates 
smaller spanwise pressure gradients on the tapered plan 
form. 



CONCLUSIONS 



1. A'swept— back wing has a lov/er total lift force and 
a higher tip loading than the corresponding unswept-back 
wing of the same plan form and aspect ratio, 

2. The changes, below the stall, caused by sweeping 
back a wing are much less pronounced for tapered plan 
forms than for constant-chord plan forms, 

3. A swept-back wing tends to stall first at the tips 
and at a lower angle of attack than the corresponding ' 
unswept-back wing. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September 16, 1941. 
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